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ABSTRACT 

We present a harmonic expansion of the observed hne-of-sight velocity field as a method to recover 
and investigate spiral structures in the nuclear regions of galaxies. Application to the emission-line 
velocity field within the circumnuclear starforming ring of NGC 1097, obtained with the GMOS-IFU 
spectrograph, reveals a three-aim. spiral in the non-circular motions, which agrees with a two-aivci 
dust spiral in the surface brightness. This nuclear spiral is consistent with a weak perturbation in 
the gravitational potential due to a two-arm logarithmic spiral, with a pitch angle of 52 ± 4° derived 
directly from the harmonic expansion of the velocity field. Next, we use a simple spiral perturbation 
model to constrain the fraction of the measured non-circular motions that is due to radial inflow. We 
combine the resulting inflow velocity with the gas density in the spiral arms, inferred from emission line 
ratios, to estimate the mass inflow rate as a function of radius, which reaches about 0.011 Mq yr^^ 
at a distance of 70 pc from the center. This corresponds to a fraction of about 4.2 x 10"'^ of the 
Eddington mass accretion rate onto the central black hole in NGC 1097, and is fully consistent with 
previous mass accretion models fitted to the observed spectral energy distribution in the nucleus 
of this LINER/ Seyfertl galaxy. We conclude that the line-of-sight velocity not only can provide a 
cleaner view of nuclear spirals than the associated dust, but that the presented method also allows 
the quantitative study of these possibly important links in fueling the centers of galaxies, including 
providing a handle on the mass inflow rate as a function of radius. 

Subject headings: galaxies: active — galaxies: kinematics and dynamics — galaxies: nuclei — galaxies: 
structure — galaxies: individual: NGC 1097 



1. INTRODUCTION 

Gas transport to the ce nters of galaxie s is still mainly 
an unsolved problem fe.g. lMartinill200^ . Since most of 
this gas is residing in a rotating disk well beyond the 
center, it is essentially a problem of angular momentum 
transport. Proposed transport mechanisms range from 
galactic interactions and bars, to nuclear bars and spi- 
rals, to stellar mass loss and disruption near the central 
black hole (BH). This range in scales also closely rep- 
resents a range in decreasing mass inflow rates, which 
in turn might be correlated with activity in the galactic 
nucleus (AGN), ranging from quasar, Seyfert, LINER to 
quiescent galaxies. However, not only is it challenging to 
(observationally) establish fueling mechanisms down to a 
few parsec from the center, time delays between changes 
in the mass inflow rate and the onset of nuclear activity 
further complicate linking them. 

It is also likely that multiple fueling mechanisms are 
important and act together. Large-scal e bars are ef- 
ficien t at transporting gas inward (e.g. lAthanassoulal 
[T99l . but the presence of an inner Lindblad resonance 
(ILR) will cause the gas to pile up in a nuclear ring, 
often clearly visible due to the intense star formation 
(see e.g. Figure [T|). The gas might continue further 
inward through nested bars (e.g. 'S hlosman et al.l Il989l : 
[Englm aicr fc Shlosman .2004), although dynamical con- 
straints on a double- barred system may prohi bit inflow 
down to the center (|Macieiewski et al] |2002[ ). More- 
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over, observational evidence of nuclear bars is scarc e 
(e.g. 'Rega n fc MulchaevI [l999l : iMartini fc Poggd I1999'). 
but should be interpreted with care since the prominent 
dust lane s in the main bar might b e absent in secondary 
bars (e.g. iShlosman fc Helleill200l) . 

On the other hand, nuclear spirals seem to be 
commonly o bserved, both in active and quiescent 
galaxies (e.g. iLaine et al.lll999i: iPogge fc M artini 2002), 
ranging f rom flocculent to g rand-design nuclear spi- 
rals (e.g. IMartini et al.l 120031 ). Whereas the former 
are suggested to for m by acoustic instabilities (e.g. 
lElmegreen et al.lll998f) . the grand-design nuclear spirals 
are thought to be the result of gas densit y waves (e.g. 
Englrnaier fc Shlosmanll2"000f l or shocks (e.g. iMaciej ewskil 
2004al lbD induced bv the non-axisymmetric gravitational 
potential of a large-scale bar. The latter two might 
be connected in the sense that the bar-driven spi- 
ral shocks trigger the gas den sity waves throughout 
the disk (|Ann fc Thakurl I2005D . which in turn seem 
to be necessary for the nucl ear spiral to be long-lived 
(jEnglmaier fc Shlosmanll2000l ). The inward extended gas 
inflow through a nuclear spiral not only depends on the 
torque of the large-scale bar, but also the gas having high 
enou gh sound speed so as to loose angular momentum 
(e.g. jEnglmaier fc Gerhardlll997l : IPatsis fc Athanassoulal 
l2000f ). as well as the presence of a central mass concen- 
tration such as a supe r-ma s sive BH to overcome a closer- 
^in (m ner ) ILR (e.g. Fukuda et al. 2000; Ann fc Thakul 

The deviations in the gas densi ty due to nuclear spirals 
are ty pically only a few per cent (jEnglmaier fc Shlosmanl 
120001 ) ■ which makes direct imaging very difficult. In- 
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stead, the obscuration due to dust thought to be asso- 
ciated with the gas overdensities is often employed, but 
when the extinction is smaU or the dust not well-mixed 
with the gas, it can lead to unclear or even missed detec- 
tions of nuclear spirals. At the same time, nuclear spirals 
induce non-circular motions in the gas, with resulting 
deviations in the observed velocity field that can be a 
significant fi'action of the underlying circular velocity. If 
the gas is ionized, the kinematics might be inferred from 
emission lines, which in general can be detected more eas- 
ily and at higher spatial resolution in the optical, than 
can be tracers of the molecular or atomic gas such as CO 
and Hi at radio wavelengths. Moreover, while imaging 
typically only yields a detection of the nuclear spiral, the 
observed non-circular motions in combination with the 
gas (over) density inferred from the simultaneously mea- 
sured (line) fluxes might be used to derive an estimate 
of the gas mass in/out-flow rate. 

An alternative approach to get a handle on the gas 
flow rates is to compute the gravitational torques from 
the observed surface brightness using a mass-to-light 
ratio conversio n calibrated against (circu lar) velocity 
measurements (jGarci'a- Burillo et all [2005h . or against 
stellar population m odels fitted to color measurements 
([Quillen et al.|[T995l ). In these cases, various assumptions 
are made, most importantly that stellar light is a clean 
(once corrected for dust obscuration) and direct (mass 
follows light) tracer of the underlying gravitational po- 
tential. Even so, this approach shows that gravitational 
torques due to non-axisymmetric structures are very ef- 
ficient at transporting gas inward while overcoming dy- 
namical barriers such as the corotation resonance. How- 
ever, when the resulting predicted mass inflow rates as 
functions of radius are compared with the non-circular 
motions in the observed gas velocity fields t he correlation 
between both is often not evident (H aan et a l.! 2009). 

The main difficulty with the non-circular motions is 
to identify the fraction that is due to pure radial flow. 
Gas on closed elliptic orbits not only contributes to non- 
circular motions (or elliptic streaming) in the azi muthal 
direc tion but also in the radial direction (Won g~et al.l 
|2004|) . However, since closed elliptic orbits are applica- 
ble strictly only to coUisionless (stellar) orbits, we expect 
that as a result of shocks and other dissipational effects 
the angular momentum of the gas will change, leading to 
in/out- flow. While the measured non-circular motions 
are often taken to be directly representative of the radial 
flow velocities fe.g. lStorchi-Bergmannll2007h . only a frac- 
tion of them are expected to be truly radial in/out-flows. 

The goal of this paper is two-fold: (i) to show that har- 
monic expansion of gas velocity fields provides a clean 
way to detect nuclear spirals and to estimate their pitch 
angles; (ii) to use simple perturbation models to con- 
strain the fraction of measured non-circular motions that 
is due to radial flow, and to estimate the correspond- 
ing mass inflow rate. In Section [2] we present the har- 
monic analysis, and in Section [3] we apply it to the ob- 
served emission-line velocity fleld within the circumnu- 
clear starforming ring of NGC 1097. We find that the 
three-arm spiral structure in the non-circular motions is 
consistent with a weak two-arm spiral perturbation of 
the gravitational potential, and use the model derived 
in Appendix |^ to constrain the inflow rate as function 
of distance from the center of NGC 1097. In Section |4] 



we discuss the corresponding spiral distortion in the sur- 
face brightness, and possible additional non-circular mo- 
tion contributions to the velocity field. Finally, we link 
the estimated mass inflow rate to the accretion onto the 
central BH of NGC 109. We summarize and draw our 
conclusions in Section [5l 

2. HARMONIC ANALYSIS 

Maps of the surface brightness (SB), line-of-sight veloc- 
ity (V), velocity dispersion (cr), and higher order velocity 
moments (often expressed in term of Gauss-Hermite mo- 
ments ft.3, hi, . . .) of nearby galaxies generally show or- 
ganized, periodic features, which can be suitably studied 
by means of a harmonic expansion. 

2.1. Harmonic expansion 

We start by dividing a map into a number of elliptic 
annuli, each with a different semi-major axis length R, 
but we assume all have the same fiattening q and share 
the same angle V'o and center (xq^i/q), so that in terms 
of Cartesian coordinates on the map (i.e., on the plane 
of the sky) 

x' = x'q + R cos ip cos ipo — 1 R sin ip sin ipo , 
v' = v'q + Rcost/jsmipo + q Rsinip cosipo- (1) 

If the map is aligned such that x' is pointing North and 
y' pointing West, then ■00 = '''/2 + F, where F is the 
common (observational) definition of the position angle 
of the major axis of the galaxy, measured from North 
through East. Next, we extract the profiles along each 
of the annuli and describe them by a finite number of 
harmonic terms n, 

n 

P = Co{R) + ^ Cra{R) COSVml) + Sm{R) SlTLTm}). (2) 
m— 1 

There are various ways to obtain the set of best-fit el- 
lipses, but given the highly non-linear nature of the op- 
timization for the best-fit parameters, this is commonly 
done in a stepwise, iterative process. We follow a s im- 
ilar approach as described in iKrainovic et all (|2006f ) — 
based on the usual ly adapted procedure in photometry 
pedrzeiewskil[l987l ) — which we briefly summarize. 

For (early-type) galaxies, the maps of SB, a and even 
higher-order velocity moments (/14, /ig, ...) each are 
to first order well approximated by a function that is 
constant along (concentric, equally flattened and orien- 
tated) ellipses, and only varies as function of the semi- 
major axis length of these ellipses. Similarly, V and odd 
higher-order velocity moments (/13, ...) each are 
well described by a similar function, but one that has an 
additional cosine variation along the ellipses. In these 
cases, the harmonic expansion in equation ([2]) truncates 
after the third-order term (n = 3). The parameters of 
the best-sampling ellipses for even velocity moments, can 
then be obtained by minimizing — I]^^]^(c^ -|- s^), 
and without ci for odd velocity moments. We first per- 
form the minimization for a grid of (fixed) q and ipo val- 
ues, where we might use external constraints such as the 
inclination i and measured position angle F together with 
an initial estimate of the center [x'^,y'^). Next, starting 
from the best-fit grid pair, we optimize for all four pa- 
rameters. 
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We sample the semi-major axis lengths as i? = i?i [fc -f- 
(1 + gY''^'^^] for k = 1, 2, 3, ... , with the initial i?i de- 
pending on the spatial resolution; for the geometric in- 
crease factor we take g = 0.1. We achieve this sampling, 
together with uniform sampling in ip, by bilinear inter- 
polation of the observed map. When fitting the observed 
line-of-sight velocity field shown in Figure [H we assume 
all motions are within the equatorial plane. We divide 
all harmonic terms by g = sinz to take into account the 
projection e ffect. For further details, including error es- 
timates, see lFathi et"all (|2005[ ). 

2.2. Spiral structure 

Once we have obtained the set of best-fitting ellipses, 
we fit the harmonic expansion of equation ^ to each of 
the corresponding profiles to obtain the harmonic terms 
Cm and Sm (up to higher order n > 3), as functions of 
R. The difference with the above expansion up to and 
including n = 3, reveals the deviations from the latter 
smooth model: for example, boxiness and diskiness in 
the SB map, and non-circular motion in the residual V 
map. 

In the case of a spiral structure as observed in 
NGC 1097, both through the dust in the SB and in the 
residual V map (|Fathi et al.l l2006f ). the deviations are 
more naturally described as a (radially varying) offset 
in the angle ■0, than through a complex combination of 
variations in the amplitudes Cm and s„i. Henceforth, 
we rewrite the harmonic expansion of equation ([2|) (in a 
mathematically equivalent expression) as 

n 

P - Ko{R) + ^"(^) ^o*^ ('"[^ - MR)]) , (3) 

m— 1 

where the amplitudes Km and phase shifts ipm are related 
to the coefficients Cm and Sm (so = 0) by 

Km=<^m+Sm ^nd tan(TO'(/'m) = — ■ (4) 

Cm 

For example, the V map of NGC 1097 in Figure [H 
reveals, after subtracting the best-fit circular motion 
Kq + Ki (R) cos ip (with Ki = ci , since we set ipi = 
si = 0), what seems to be a three-arm spiral structure, 
K2{R) cos{3['4' — '4'3{R)]). Here, ipsiR) traces the spiral 
arms as function of radius R. 

2.3. Weakly perturbed gravitational potential 

An axisymmetric gravitational potential in a frame 
that co-rotates with a weak perturbation of harmonic 
number m can be written in polar coordinates as 

$(i?, 0) = $o(i?) + ^m{R) cos (m [0 - (j)miR)]) , (5) 

where 4'm is the phase of the perturbation. Through 
Poisson's equation it follows that the corresponding sur- 
face mass density exhibits a harmonic m distortion. 

To derive the li ne-of- sight velocity, we follow 
ISchoenmakers et aD (|1997t ). As described in Ap- 
pendix[Al we extend their collisionless analysis by includ- 
ing radial damping in the equations of motion to take 
into account the dissipative nature of gas. We assume 
that the gas moves on closed loop orbits in the equa- 
torial plane, which we observe at an inclination i away 
from its normal and at an (aziniuthal) angle 4>ios- Given 



a point (i?, "0) in the equatorial plane, the projection of 
the azimuthal and radial velocity onto the corresponding 
line-of-sight yields 

Vios = sinz [v^{R,ip) cosip + vr{R,iI}) sinV'] , (6) 

where = cf) — (^in^ + 7r/2 is zero on the h ne of nodes (see 
also Fig. 1 of ISchoenmakers et al.|[l997f ). To first order 
the solutions of the equations of motion in the perturbed 
gravitational potential of equation ^ yield 

vr{R,^)=Vc{R) [cr cos rmp + SRsinmip] , (7) 
Vci,{R, ip) — Vc(R) [1 -I- C0 cos mip + sinm-0] , (8) 

where v1 = _Rd<i>o/di? is the circular velocity and cr^ 
sr, C(j, and are functions of R given in Appendix lAl 
Substituting these solutions into equation ([S]), we obtain 

Vios = K cos tp 

+ Cm-1 cos(m — l)ip + Sm-1 sin(m — l)xp 

+ Cm+i cos(to + 1)iIj + Sm+1 sin(m -1- l)i/;, (9) 

with Cm±i = Kt(c0 T Sfl)/2 and Sm±i = 14(s0 ± c_r)/2, 
and V* = Vc{R)sini the circular velocity in projection. 
We th us find, as concluded before by Schoenmakers et all 
(1997*) and already qualitatively inferred by ICanzianI 
(1993) that, if the gravitational potential has a pertur- 
bation of harmonic number m, the line-of-sight velocity 
field contains an m — 1 and an m -f 1 harmonic term. 

2.4. Pitch angle 

How loosely or tightly wound a spiral is can be quan- 
tified via its pitch angle, C, which, at a given radius, 
R, measures the angle between the tangent of the spi- 
ral arm and a circle with radius R in the plane of the 
disk. Inverting the phase-shift, (f>{R)^ of a spiral, we can 
parameterize an arm of the spiral as 

X — R{(j)) COS0 and y = R{(f>) sin0. (10) 

One can then show that the pitch angle C is given by 



which is positive (negative) if the spiral curves anti- 
clockwise (clockwise) with increasing radius. The smaller 
the pitch angle, the more tightly the spiral is wound, with 
C = a circle, while ( = ±7r/2 corresponds to a straight 
line. A specific case that is often encountered in nature 
is that of a logarithmic spiral 

1 R 

(t>{R) = — In — ^ R{(t)) = an exp(6o 0), (12) 
bo oq 

with constants oq and bo- From equation (|lip . we find 
the well-known property that the logarithmic spiral has 
a constant pitch angle ^ — tan^^ bo. 

In the case of an (m ± l)-spiral in the line-of-sight ve- 
locity field, the corresponding pitch angle, Cm±i7 follows 
directly from the phase-shift, ipm±iiR), in the harmonic 
expansion in equation The pitch angles Cm-i and 
Cm+i generally take different values and bracket the pitch 
angle C„j of the m-spiral perturbation in the gravitational 
potential that caused them. To show this, we start from 
equations (|A12p and (jA13p in Appendix |^ and rewrite 
the coefficients Cm±i and Sm±i as 

Cm±l^Km±l COsimipm - 6m±l), 
Sm±l=Km±lSin{mipm- Om±l), (13) 
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so that after substitution into equation Q we obtain 

mifjn- 0„i±l ^ {m±l)iprn±i, (14) 

where ipm — <i>m — 4>\o& + 7r/2. This hnks the phase-shift 
(j)miR) of the m-spiral perturbation in the gravitational 
potential with the phase-shifts ipm±i{R) of the (m ± 1)- 
spirals in the linc-of-sight velocity field. The correspond- 
ing pitch angles are related as 

m cot(„i = (m ± 1) cotCmii - d9,n±i/d\iiR. (15) 

In general, Km±i and 0m±i depend in a rather complex 
way on the gravitational potential, but as we show in 
Appendix [Al thev possess some generic properties. 

First, the amplitude Km+i is larger (smaller) than the 
amplitude Km-i outside (inside) the corotation radius 
i?CR) and equal to it at ii!cR- This implies a transition 
in the line-of-sight velocity field at Rqki going from the 
(m + l)-spiral dominating outside i?cR to the (m — 1)- 
spiral dominating inside Rgk (see also Fi gure E I . This 
is also concluded by Schocnmakers et aL (Il997[l for the 
collisionless case and earlier by ICanzianI (|1993f ) for the 
less general case of a tightly wound spiral in the lin ear 
density- wave theory (see also lCanzian fc Allenlll997l ). 

Second, 9m±i typically varies much less with radius 
than the spiral phase-shifts, so that d6'm±i/dlni? in 
equation (|15p is relatively small. As a result, we can 
estimate the pitch angle Q„i of the m-spiral perturbation 
in the gravitational potential from the pitch angles Cm-i 
and/or Cm-(-i of the (m— l)-spiral and (TO-l-l)-spiral in the 
observed line-of-sight velocity field, without constructing 
a full dynamical model (see also the last two panels of 
Figure [5]) . 

2.5. Radial flow velocity 

Nuclear spirals in principle provide a mechanism to 
transport gas from kpc scales, where it often stalls inside 
a nuclear ring, into the center of the galaxy. Still, as 
mentioned in Section [1] if a nuclear spiral (or a nuclear 
bar) is due to gas moving on closed elliptic orbits this 
results in non-circular motions, also referred to as ellip- 
tic streaming, but not necessarily in net inflow towards 
and/or outflow away from the center. However, unlike 
stars, gas is not collisionless and its orbits interact and 
ex- c hange angula r momentum leading to net radial flows 
(e.g. lWadalfl99l . 

Henceforth, in the analytic models in Appendix El 
we assume a weak perturbation in the gravitational po- 
tential causing gas to deviate from circular onto ellip- 
tic orbits, while taking into account its dissipative na- 
ture via radial damping. These analytic spiral models 
are an extension of the analytic bar models introduced 
by (jWadal [l994j) . who showed that they describe well 
the gas behavior seen in hydrodynamical simulations. 
The radial damping causes the gas to lose/gain angular 
momentum inside/outside the corotation radius, which 
nicely matches the angular momentum transfer due to 
the torque from the bar potential in numerical simula- 
tions. The amount of radial damping, controlled through 
the dimensionless parameter'^. A, thus provides a handle 
on the amount of net radial flow that is needed to explain 

^ Note that o ur A is the same as the dimensionless parameter A 
in lWadal iITQQI '). while we define A = 2XKm(Q. -Q.p). 



the observed non-circular motions in addition to elliptic 
streaming. 

In this way, a rather straightforward estimate of the net 
radial flow velocity, Ufiow , can be obtained by comparing 
the radial velocity vr (equation lA5|l of the analytic model 
that includes radial damping (A > 0) with the analytic 
model without radial damping (A = 0). At a given radius 
i?, the maximum radial velocity is given by 

VR^^a.yi^rn{n-Vlp)Ri-^--—^\ . (16) 

Here, A and B only depend on the weak perturbation as 
given in equation (|A2p , while 

= -m^iVL-VLpf, A^2\Km{fl-np), (17) 

are functions of the angular frequency 0(i?) and the 
epicycle frequency k{R) of the axisymmetric gravita- 
tional potential, as well as the harmonic number m and 
pattern speed flp of the weak perturbation. Since with- 
out radial damping (A = 0) all radial motion is due to el- 
liptic streaming, we subscribe a fraction | A|/(A^-|-A^)^/^ 
of the radial velocity to elliptic streaming, leaving as an 
estimate of the radial flow velocity 

/ A2 

"flow = [ ^2^f^2 J '"R- (18) 

Even though the analytic models neglect possible non- 
linear effects, they can capture most features of observed 
non-circular motions (as we show next in the case of 
NGC1097), and at the same time provide an estimate 
of the fraction of the observed non-circular motions that 
is due to net radial flow in addition to elliptic stream- 
ing. Note that radial flow here does not mean that the 
gas is following pure radial orbits with zero angular mo- 
mentum, whi ch would contribu te only to the harmonic 
term si fe.g. iWong et al.|[2004D . Instead, the gas is ex- 
pected to gradually spiral inward/outward as it has both 
azimuthal and radial velocity components. In case the 
angular momentum loss/gain is driven by a weak grav- 
itational potential perturbation with harmonic number 
m, this results in a contribution to Cm±i and Sm±i- 

3. NUCLEAR SPIRAL IN NGC 1097 

We apply the above harmonic analysis to the observed 
emission-line velocity field within the circumnuclear star- 
forming ring of NGC 1097. We recover in the non-circular 
motions a spiral structure and infer its pitch angle di- 
rectly from the harmonic components. Next, we use a 
spiral perturbation model to constrain the radial inflow 
velocity and combine this with the gas density in the nu- 
clear spiral to estimate the mass inflow rate as function 
of distance from the center of NGC 1097. 

3.1. Non-circular motions 

NGC 1097 (ESQ 325-58) is a nearby (distance 
14.5 Mpc, so 1" ~ 70 pc) LINER/Seyfert 1 host with a 
strong, ~ 16 kpc long, bar and a ~ 0.7 kpc in radius 
circumnuclear starforming ring (Figure [1]). In this Sb 
galaxy dust can be traced within the large-scale spiral 
arms out to a (outer Lindblad resonance, OLR) radius 
of ~ 14 kpc, in prominent lanes along the bar, and con- 
tinuing within the nuclear ring as a spiral structure down 
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Fig. 1. — Left: VLT/VIMOS color composite image of NGC 1097, showing the large-scale spiral arms and bar with prominent dust lanes 
reaching down to the circumnuclear ring (credit: European Southern Observatory). Middle: VLT/NACO adaptive optics color-composite 
image of the circumnuclear ring region (credit: European Southern Observatory), with the footprints of the observations with the GMOS- 
IFU spectrograph. Right: HST/ACS structure map of the same circumnuclear region with the wavelet map of lLou et al.l ((2001) overplotted 
in color, with increasing intensity from blue to red. The solid magenta curves show the two-arm nuclear spiral with pitch angle 52°, which 
we predict based on the three-arm spiral structure in the velocity field. The dashed magenta curves indicate the two additional arms in 
case a m = 4 spiral perturbation with the same pitch angle would be present (see Section |4. II for further details). 
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Fig. 2. — The first panel shows the line-of-sight [Nil] emission-line velocity field within the circumnuclear ring in NG C 10 97, obtained 
with the GMOS-IFU spectrograph. The other panels show the results of applying the harmonic expansion of Section 12.11 the best-fit 
circular motion, the remaining non-circular motion, and the harmonic reconstruction from the sum of the latter two. At the right-hand 
side of each map, the (linear) velocity scale is indicated by the color bar, and the limits are given below each bar. The black cross indicates 
the (kinematic) cent er. T he overplotted logarithmic spiral with pitch angle of 52° traces the three-arm spiral structure in the non-circular 
motions (see Section |3. II for further details). 
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Fig. 3. — Harmonic expansion of the observed line-of-sight velocity field of NGC 1097. After fixing the position angle F and systemic 
velocity Vaya to their mean values indicated by the dotted line in, respectively, the left-top and middle-top panels, the harmonic coefficients 
Cm and Sm are extracted as functions of radius R (in arcsec). The coefficients have been divided by g = sin j to correct for the inclination 
i = 35°. The dotted vertical line in each panel indicates the radius R = 1", within which the measurements are less accurate due to the 
presence of dust. In the bottom panels, the phase-shifts ipm defined in equation lj4)l are plotted against I n R , so that a spiral structure shows 
up as a smooth variation, which becomes linear in case of a logarithmic spiral as shown in equation I I12I I. In this way, the bottom-right 
panel shows that the structure in the non-circular motions is consistent with a three-arm logarithmic spiral with indicated pitch angle 1^3 
from the slope of the thick solid linear fit. The indicated estimate of the pitch angle ^2 of the inferred two-arm spiral perturbation in the 
gravitational potential follows from equation I I15I I , if we neglect the radial variation in the additional term 63 . The latter is indeed expected 
to be small as indicated by the dot-dashed curve, which shows 63 based on an analytical two-arm spiral perturbation model. Taking into 
account this term yields the model prediction for ^3 indicated by the dashed curve, which is nearly indistinguishable from the fitted thick 
solid line. Similarly, in the bottom-left panel, the dot-dashed and dashed curves show the model prediction of 81 and ■ipi , respectively. In 
line with the measurements, ipi is small because the ci coefficient in tanipi = si/ci includes the dominating circular velocity Vc- The latter 
is shown as a black dashed curve in the top-right panel, whereas including the non-circular motions due to the modeled two-arm spiral 
perturbation, yields ci indicated by the blue dashed curve. Similarly, the model prediction for si is given by the red dashed curve in the 
middle- left panel, and those for C3 and C3 by the blue and red dashed curves in the middle-right panel. (See Section l3.2l and Appendix lAl 
for further details on the spiral perturbation model.) 
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to < 3.5pc from the center (|Lou et al.ll2001l : lPrieto et all 
|2005| ) . Non-circu lar motions associ ated with this nuclear 
spiral structure (|Fathi et al.1 12006[ ). indicate a possible 
mechanism to drive gas from kpc scales down to a few pc 
from the ce nter, where a double-peak ed broad Ha emis- 
sion profile (jStorchi-Bergmann et al.lll993,) indicates the 
p resence of a s uper- m assive BH. 

iFathi et all (|2006[ ) describe in detail the observa- 
tions and reduction of the two-dimensional spectroscopy 
of NGC 1097 obtained with the GMOS-IFU on the 
Gemini South Telescope (GS-2004B-Q-25, PI: Storchi- 
Bergmann). Three pointings within the nuclear ring re- 
gion (Figure [1]) provided 1500 individual spectra cover- 
ing 5600-7000 A at a velocity resolution of 85 km s ^ and 
with a spatial sampling of O'.'l. 

In the first panel of Figure [21 we present the [Nil] 
emission-line velocity field covering the inner 0.5 x 
l.Okpc. A systemic velocity of fgys = 1188kms~^ 
(equivalent to a constant cq term) has been subtracted. 
For an inclinati on i = 35° (flattening q = 0.82; 
IFathi et al.ll2006[ ). and position angle of F = 141°, the 
other panels show the results of applying the harmonic 
expansion of Section 12.11 From left to right: the best- 
fit circular motion V^irc — ci cos ip, the remaining non- 
circular motion, and the harmonic reconstruction from 
the sum of the latter two. 

In Figure [H the first two panels show the adopted po- 
sition angle and systemic velocity. The next four pan- 
els show the coefficients Cm and Sm (in kms~^) for the 
first three harmonic terms, as function of the (depro- 
jected) radius R. The bottom three panels, show the 
phase shifts ipm (in degrees) for the alternative formula- 
tion of the harmonic expansion given in equation ((3|) , as 
functions of the (natural) logarithm of R. 

The apparent three-arm spiral structure in the non- 
circular motion (third panel in Figure [2]) is confirmed 
by the significant amplitudes of C3 and S3, and in par- 
ticular by the the smooth variation of the corresponding 
phase-shift ips with radius. Excluding the uncertain mea- 
surements within R = 1" (vertical dotted line) and the 
single measurement at the edge of the map, the remain- 
ing measurements indicated by solid circles show a clear 
linear relation between Sips and In R. This is consistent 
with a logarithmic spiral defined in equation (I12p . with 
fitted ao = 1.82 ± 0.33" and 60 = -1.95 ± 0.27. The 
slope provides a direct and robust measurement of the 
pitch angle of ^3 — 63 ± 3°. This logarithmic spiral, 
overplotted in the third panel of Figure [21 indeed closely 
traces the three-arm spiral structure in the non-circular 
motions; possible reasons for the deviations are discussed 
below in Section [H 

The positive sign of the pitch angle indicates the nu- 
clear spiral is curved anti-clockwise with increasing ra- 
dius, equivalent to the spiral arms extending outwards of 
the large-scale bar, as can be seen in the left panel of Fig- 
ure[T] Since the global rotation in NGC 1097 is clockwise 
from the observed velocity field, and the North-East is 
the "far side" from being more obscured, it follows that 
both the large-scale and nuclear spirals are trailing. 

3.2. Two-arm spiral perturbation 

The above three-arm spiral structure in the non- 
circular motions is consistent with a perturbation in the 
gravitational potential due to a m = 2 harmonic spi- 



ral (see also Section 12. 3[) . From equation (fT5)) it then 
follows that the pitch angle ^2 of this two-arm spiral per- 
turbation in the gravitational potential follows from the 
measured pitch angle Cs of the three-arm spiral in the 
velocity field as 2 cot C2 = 3 cot (3 — dO^ /d\nR. We show 
below that, as expected, the later term is small, so that 
for a measured Cs — 62 ± 3°, we calculate C2 — 52 ± 4°. 

As illustrated in Figure[5]of AppendixIXl the true value 
of C2 might differ slightly depending on the details of 
the gravitational potential perturbation. We construct 
a model for the nuclear spiral in NGC 1097 to get a 
handle on the latter difference as well as to constrain 
the fraction of the observed non-circular motions that is 
due to net radial inflow in addition to elliptical stream- 
ing. We use the analytic solutions of Appendix [^ for 
gaseous orbits in an axisymmetric gravitational poten- 
tial, $0(^)1 that is weakly perturbed by a logarithmic 
m = 2 spiral with pitch angle (2- Note that these ana- 
lytic spiral models are based on linearized equations of 
motion under the epicycle approximati on, but are not 
restricted to tightly w ound spirals (e.g. iLin et al.|[l969l : 
ICanzian fc Allenlll997D . In this way, we show below that 
the observed loosely wound nuclear spiral in NGC 1097 is 
still consistent with a density wa ve, and not necessa rily 
driven by shocks as suggested by'Davies et al.' (20 091). 

We adopt the power-law model (Evans & de Zeeuwl 
Il994[ ) with axisymmetric gravitational potential 
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so that Vc{Rc) = fo/2 at the core radius Rc- The param- 
eter (3 controls the logarithmic gradient of the rotation 
curve at large radii: (3 < Q rising, (3 = Q flat, and (3 > 
falling. The three gravitational potential parameters vq , 
Rc and (3 are set by comparing the corresponding circular 
velocity Vc{R) with the measured radial profile of ci, tak- 
ing into account the non-circular contribution due to the 
m = 2 spiral perturbation. In the third panel of Figure[3l 
the black dashed curve shows Vc{R) for vq — 275kms~^, 
Rc = 1.2", and (3 = —0.6. Except for the three values 
within R = 1", this simple power-law model is a good 
representation of the ci measurements, once the non- 
circular contribution derived below is added, as indicated 
by the blue dashed curve. The advantage of such a sim- 
ple analytic representation of the gravitational potential 
is that it makes all subsequent calculations concerning 
the perturbation very convenient. 

For the amplitude of the gravitational potential per- 
turbation we assume ^2{R) = £p^o{R), with constant 
strength e^, while the phase-shift is given by 4'2{R) — 
cot C2 ln(i?/ao). Here, oq = 1.8" from the above fit to 
the m' = 3 harmonic terms in the line-of-sight velocity 
field, while the corresponding approximation (^2 — 52°, 
is taken as the initial value for the pitch angle. The 
three additional free parameters are the pattern speed. 
Sip, of the spiral perturbation, the azimuthal viewing 
angle, ^los, and finally the amount of radial damping. 
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A. Equations (|A12[) and (jA13|l then provide predictions 
for the non-circular motion contribution in terms of the 
harmonic coefficients ci, si, C3, and S3, which we com- 
pare with the corresponding measured radial profiles for 
NGC1097 in Figure [3 

The flattening of the Hi rotation c urve in the outer 
parts of NGC 1097 (jSofue et all I1999D implies a nearly 
constant circular velocity Vc — 300kms~^. With the 
angular and epicycle frequencies approximately given 
by r2 ~ ~ Vc/R, the OLR being at i?oLR — 

14kpc yields a pattern speed flp ~ il. + k/2 ~ 
35kms^^ kpc^^. The latter places the corotation ra- 
dius, RcR — Vc/^p — 8.6 kpc, or ~ 10% beyond 
the extent of the large-scale bar, consis tent with nu- 
merical simulations of "fast bars " (e.g. lAthanassoulal 
Il992t iDebattista fc Sellwoodl[2000l ). and me asured pat- 
tern speeds in simil a r galaxies (e.g. lA gucrri et al.l 120031 : 
iGerssen et al.l l2003t iRautiainen et al.i i2008i) . Assum- 
ing that the two-arm nuclear spiral as a gas den- 
sity wave is being driven by the large-scale bar 
(|Englmaier fc Shlosman|[2000f ). we adopt the same value 
for the pattern speed of the perturbation. Since we con- 
sider the harmonic coefficients well within the corotation 
radius, they are not sensitive to, and hence do not con- 
strain Qp. 

On the contrary, a small change in the azimuthal 
viewing angle already causes a significant radial offset 
in the harmonic coefficients, so that we need 0ios — 
(1.05 ± 0.05) 7r/2. Next, too little radial damping re- 
sults in too small amplitudes of the coefficients ca and 
S3 with respect to si and the non-circular contribution 
to ci. Specifically, since si and S3 are of similar ampli- 
tude and shape over most of the radial range, significant 
radial damping with A > 1 is needed. As described in 
Section I^TSl and discussed below in Section 1^751 this im- 
plies that net radial flow makes up most of the intrinsic 
radial velocity. Matching the amplitudes of the harmon- 
ics terms yields a strength of the gravitational potential 
perturbation of Cp ~ 0.15. Finally, values for the pitch 
angle of the spiral perturbation that are in the range of 
C2 — 52 ± 4° , approximated above from the pitch angle 
(3 of the three-arm spiral in the velocity field, yield pre- 
dictions for the harmonic coefficients that are consistent 
with the measured harmonic coefficients. The effect of 
the additional term d6'3/dlni? is indeed small, so that 
2 cot C2 ~ 3 cot C3 provides a robust measurement of C2 • 

This is further illustrated in Figure [3l where in ad- 
dition to the measured harmonic coefficients, we show 
with dashed curves the predictions of the above two-arm 
spiral perturbation model with pitch angle C2 = 52° , az- 
imuthal viewing angle (f>\os — 94.5°, and radial damping 
parameter A = 2. Both 9i and 62, shown as dot-dashed 
curves in respectively the bottom-left and bottom-right 
panel, indeed vary only little with radius. Similarly, the 
dashed curve in the bottom-right panel shows the pre- 
dicted relation between SV's and Ini? while taking into 
account ^3, which is nearly indistinguishable from the fit- 
ted solid line. The predicted harmonic terms match well 
the C3 and S3 measurements in the middle-right panel, 
and rather well the si measurements in the middle-left 
panel, except for those at larger radii. However, the lat- 
ter measurements are quite uncertain, and might well be 
disturbed by the nuclear ring. The blue dashed curve 



in the top-right panel is the spiral model prediction for 
ci, which apart from the innermost measurements, nicely 
traces the measured rotation curve. 

Due to the significant non-circular motion contribu- 
tion the latter is different from the circular velocity 
of the (power-law) axisymmetric gravitational potential 
shown as the black dashed curve. The corresponding 
r2(i?) — k{R)/2 black dot-dashed curve is still well above 
the pattern speed f2p = 35kms~^kpc~^ (solid horizon- 
tal line) at radius ^ 8", which places the ILR well be- 
yond the nuclear ring radius of ~ 0.7 kpc. If indeed 
the radius o f nuclear rings is set by the location of 
the ILR (e.g. iButa fc Combe£lll996l ). this indicates that 
the nuclear ring in NGC 1097 has migrated inward. In- 
ward migration has also been suggested for the nuclear 
ring in NGC 4314 (Be nedict et al.ll2002D. has been seen 
in hydrodynamic simu lations (e.g. iFukuda et all 120001 : 
iRegan fc Teubenll2003f) . and might be a consequence of 
shep herding of the gas ring by st ar clusters that formed 
in it (|van de Ven fc Changj 120081) . 

3.3. Mass inflow rate 

We found that a spiral perturbation model with only 
elliptic streaming, in which gas is moving on closed ellip- 
tic orbits without radial damping, cannot explain the ob- 
served non-circular motions in NGC 1097. In Section [^31 
we showed that the amount of radial damping required 
provides an estimate of the radial flow velocity vaow in 
terms of the radial velocity vr of the spiral model which 
matches the observed harmonic terms. Since in the case 
of NGC 1097 a large amount of radial damping A > 1 
is needed, the ratio Wflow/wfl — jAj/VA^ + A^ in equa- 
tion p8)) is close to unity. In other words, the contribu- 
tion to the radial velocity in the nuclear spiral is predom- 
inantly due to net radial flow, which is directed inward to 
the center of NGC 1097, as we are well within the coro- 
tation radius. For the same spiral model as shown in 
Figure [3] with the dashed curves, the corresponding ra- 
dial inflow velocity ^inflow is shown in the top-left panel 
of Figure [4] with a solid curve. Even though nearly iden- 
tical to the radial velocity vn, it still is only a fraction of 
the total non-circular motion (dashed curve), which also 
includes the azimuthal velocity. 

To turn this radial inflow velocity into a constraint on 
the mass inflow rate, we furthermore need to know the 
gas density in the nuclear spiral, as well as the geometry 
or the area through which the gas is flowing. The flux 
ratio of the [Sii]AA6716.4,6730.8 emission-line doublet in- 
cluded in the spectral range of the GMOS-IFU observa- 
tions allows us to constrain the mean electron density 
Tie over the observed field. Taking all values together 
yields an average flux r atio of 1.003ib0.031, wh ich, adopt- 
ing the prescription bv lShaw fc Dufoud (IT9 9I and using 
the atomic p aram eters compiled by lMendo za (1983) and 
IOsterbro"c3 ()1989l ). corresponds to He = 600 ± 77 cm"'^ 
assuming a mean electron temperature of lO'^K. Even 
though the flux ratio map might indicate a spiral struc- 
ture similar to that in the non-circular motions, the 
signal-to-noise is not high enough to quantify the gas 
(over)density in the spiral arms from it. Instead, we 
apply the same elliptic annuli used to extract the har- 
monic coefficients to the flux ratio map. From the distri- 
bution within each annulus, we derive Ue as a function 
of radius, resulting in the dashed curve in the top-right 
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Fig. 4. — Estimate of the mass inflow rate as function of distance R (in pc) from the center or NGC 1097. Top-left: The solid curve shows 
the radial inflow velocity inferred from a spiral perturbation model matched to the harmonic expansion of the velocity field. The radial 
inflow velocity is only a fraction of the non-circular motions indicated by the dashed curve. Top-right: The solid circles show the electron 
density Arte in the nuclear spiral arms estimated from the variation in the flux ratio of the [Sll] AA6716.4,6730.8 emission-line doublet. This 
(over)density is typically a factor ten smaller than the average electron density ne as indicated by the dashed line. Bottom-left: The solid 
circles show the adopted scale height h, assuming a marginally stable disk resulting in the dashed curve, but with the relative width R/ (2m) 
of the m = 2 spiral mod el in the equatorial plane as an upper limit. Bottom-right: Combining the solid-curve values of the first three 
panels into equation l|2H l results in the shown mass infiow rate (in MQyr~^). The upper dashed horizontal line indicates the Eddington 
accretion rate onto the central black hole in NGC 1097 with a mass Mbh = 1.2 X 10* based on the central stellar velocity dispersion. 
The lower dotted horizontal line at a factor 0.01 of the Eddington accretion rate indicates the approximate transition from a Seyfertl to a 
LIN ER active galactic nucleus, and is thought to be the rate below which mass accretion becomes radiatively inefficient. (See Sections 13.31 
and 14.31 for further details.) 



panel of Figure[3] (divided by 10 for illustrative purposes). 
Next, we assume that the width of the distribution, in- 
dicated by the solid curve, is driven by the density wave 
contrast, and hence yields an estimate of the electron 
(over)density Arte within the spiral arms. Assuming a 
50% lower (higher) electron temperate of 0.5(1.5) x 10'* K, 
leads to an increase (decrease) in the electron densities 
by a factor of about 20%. In all cases, the density con- 
trast is of the order of 10 %, in agreeme nt with hydro- 
dynam ic simulations (e.g. lEnglmaier fc S hlosman 2000; 
Macieie wski 2004b) and if-band imaging presented by 
Davies et al.i (|2009(1 . Using the proton mass and a factor 
1.36 to account for the presence of Helium, we convert 
Arte to a gas mass overdensity Apgas- 



For the geometry we assume in agreement with the spi- 
ral perturbation model that both the radial flow velocity 
and gas mass overdensity, at each radius i? in a disk with 
scale height h, vary as a sinusoidal function in azimuthal 
angle "0. Integrating over the ip values for which the ra- 
dial flow velocity is positive, the mass inflow rate then 
reduces to 

M = m Winflow Apgas ttR^ 4 ' (21) 

R 4m 

for gas of density Apgas flowing in the m spiral arms 
towards the center at a velocity Winflow To estimate 
the scale height h, we sta rt from a marginally stable 
disk with Toomre's (|1964f ) Q ~ 1, and substitute in 
Q — CsK/TrGSgas a constant sound speed Cs — lOkms"*, 
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the epicycle frequency n of the power-law axisymmetric 
gravitational potential, and Sgas — Pgash, with pgas in- 
ferred from the average electron density rie- The dashed 
curve in the bottom-left panel of Figure [4] shows h/R 
as function of radius. Clearly, towards the center this 
leads to an unrealistically high scale height, so that we 
constrain h/R to be not larger than the relative width 
l/(2m) of a spiral arm, resulting in the solid curve. Fi- 
nally, in the bottom-right panel of Figure 01 we present 
the mass inflow rate M (in M0 yr"-'^) as function of ra- 
dius. 

4. DISCUSSION 

We have argued for a two-arm spiral perturbation in 
the gravitational potential as the source of the three-arm 
spiral structure in the velocity field of NGC 1097. We 
verify the implied two-arm spiral distortion in the surface 
brightness, and discuss possible causes for contributions 
from additional harmonic terms. Finally, we link the 
mass inflow rate to the accretion onto the central BH. 

4.1. Distortion surface brightness 

An m — 2 spiral perturbation in the gravitational po- 
tential that would explain the three-arm spiral structure 
in the velocity field of NGC 1097, implies an m = 2 
spiral distortion in the surface mass density and hence 
a two-arm spiral structure in the surface brightness. 
Even if the latter distortion is too weak to measure di- 
rectly, it might show itself through obscuration by cor- 
related dust. Indeed, structure maps of N GC 1097 (e.g. 
iPogge fc Martinill2002t iMartini et gd] 120031) show spiral- 
like features, but they ar e not evidently a two-arm nu- 
clear spiral. iPrieto et al.l ()2005l ) note in high-resolution 
VLT/NACO infrared images a central spiral with a three- 
arm symmetry, although one of the three arms does 
not seem to continue towards the nuclear ring, but in- 
stead splits into a number o f spiral filaments. On the 
other hand, iLou et all (|2001h use wavelet analysis on a 
HST/WFPC2 image to reveal a clear two-arm spiral struc- 
ture, but additional spiral features cannot be ruled out. 

In the ri ght panel of Figu re [H we reproduce the wavelet 
map from (iLou et al]|200lL t heir F ig. 1) on top of a struc- 
ture map from iFathi et al.l ()2006l ). The latter structure 
map is based on the Ric hardson-Lucy image restoration 
technique (.Snvder et al . 1993), using a multi-step convo- 
lution of the HST/ACS high resolution camera FR656N 
image with a tw o-dimensional PSF m odel constructed 
using Tiny Tim (jKrist fc Hookl 119971 ). The solid ma- 
genta curves show the two-arm spiral with pitch angle 
^2 = 52°, which we predict based on the three-arm spiral 
in the velocity field (Section I3.2|l . The two open arms 
trace well the spiral structures in both the wavelet and 
structure map, except closer to the nuclear ring when, at 
least in the northern part, a more tightly wound spiral 
seems needed. 

The dashed magenta curves indicate the two addi- 
tional arms in case of an m = 4 spiral perturbation with 
the same pitch angle. Such higher-order even harmonic 
terms can result from non-linear coupling of modes, and 
— although smaller in amplitude then the m = 2 spiral 
perturbation — might give rise to possibly additional spi- 
ral features. Furthermore, a spiral perturbation driven 
by the large-scale bar is not necessarily restricted to an 



m = 2 harmonic term'*. Still, as long as the perturbation 
is bi-symmetric all resulting harmonic terms are even. 

Only in case of an asymmetric perturbation, or lopsid- 
edness, are contributions of m = 1 or higher-order odd 
harmonic terms to the surface brightness possible. Since 
the stars make up nearly all the mass in the center of 
NGC 1097, such a deviation from bi-sy mmetry, e.g. du e 
to an eccentric nuclear disk as in M31 (|Tremainelll995l) . 
or a dwarf elliptical companion, should reveal itself in the 
stellar kinematics. However, the central stellar kinemat- 
ics of NGC 1097 do not show such deviations nor any (spi- 
ral) s tructure in the non-circular motions (jPavies et al.l 
|2009( ). in support of the n uclear spiral b eing a pertur- 
bation in the gas. Instead, iPrieto et al.l (|2005[ ) propose 
that a difference in strength between the two principal 
shocks in the large-scale bar give rise to the lopsided- 
ness. This only works in their scenario in which the 
shocks, and not the bar itself, create the nuclear spiral. 
However, while the shocks might be the trigger, it seems 
that density waves in the gas, driven by the bi-symmetric 
bar, are necessary for the nuclear spiral to be long-lived 
(|Englmaier fc Shlosmanll2000t lAnn fc Tha kur 2005). Al- 
ternatively, the typically asymmetric (projected) dust 
distribution might mimic odd harmonic terms in the sur- 
face brightness, and possible even give rise to even har- 
monic terms in the velocity field, as we discuss next. 

4.2. Additional non-circular motion 

The coefficients Cm and Sm, that provide the phase- 
shift ipmiR): may be "contaminated" by additional con- 
tributions to the non-circular motion, but in general they 
are not expected to result in a smooth variation of the 
phase-shift with radius, as seen for a spiral structure. 
An exception is the m = 1 harmonic term, since the 
measured coefficient ci also incorporates the circular ve- 
locity as UcSini, which in general dominates over the 
non-circular motions c"*^ and si. This results in a phase- 
shift -01 that is everywhere close to zero, and due to the 
degeneracy in ci there is little hope of constraining a 
two-arm spiral perturbation from its contribution to the 
m — I harmonic term in the velocity field. However, 
we can predict its pitch angle Ci and phase-shift 
from the (logarithmic) spiral inferred from the contri- 
bution to the TO = 3 harmonic term: tan Ci — 3 tan C3 
and tp^'^{R) ~ 3ip^'^{R). As a result, we might use 
cf^ = si cot tpi'^ as an estimate of the non-circular con- 
tribution to ci. The effect is shown in the top-right 
panel of Figure [31 where the black dashed curve is Vc 
of the (power-law) axisymmetric potential, which after 
taking into account c5"^ yields the blue dashed curve that 
matches the measured ci coefficients indicated by the 
blue diamonds. This provides a novel way to correct for 
non-circular motions, which otherwise might, for exam- 
ple, lead to an underestir nation of the (inner) slop e of 
the mass distribution fe.g. iHavashi fc Navarroll2006l ) . 

Figure [3] shows that, besides m = 1 and to = 3 
harmonic terms in the velocity field expected from a 
m = 2 perturbation in the gravitational potential, the 
non-circular motions also seem to contain an m = 2 har- 

* For example, the axisymmetric power-law potential in equa- 
tion l|19| l. being perturbed by replacing the radius R by R = 
-I- {y/q)^ with (J < 1, creates besides m = 2 also higher-order 
even harmonic terms. 
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monic term. There are several effects that might (partly) 
cause this addit ional contribution. 

As shown by ISchoenmakers et al] (|1997l their equa- 
tion 7), an error in the (kinematic) center results in a 
spurious m — and m — 2 contribution to the line-of- 
sight velocity as 

We see that the effect on C2 and S2 vanishes if a = 1, 
i.e., when the circular velocity curve increases linearly 
with radius, Vc oc R. As in most galaxies, the latter 
is also the case in the inner region of NGC 1097, but 
still C2 is significantly positive within the central < 2". 
Not surprisingly, varying the kinematic center (in the 
process of finding the best-fit set of ellipses as described 
in Section 12. ip does not remove the m = 2 harmonic 
contribution, but in contrast makes the overall fit worse. 
This makes a significant effect due to an incorrect center 
unlikely. 

Based on the analysis in Section 12.31 we expect a sim- 
ilar contribution of both m — Q and m ~ 2 harmonic 
terms to the velocity field from an m = 1 distortion of 
an axisymmetric distribution. As mentioned above in 
Section [521 it is unlikely that the gravitational potential 
itself is lopsided since the stars that dominate in mass 
do not show any such signature. Also, the nuclear spi- 
ral itself is expected to be bi-symmetric if it is indeed a 
gas density wave driven by the large-scale bar. Still, the 
phase-shift 2ip2 in the bottom-middle panel of Figure [3] 
seems to vary smoothly, and even close to linearly as a 
function of Ini? within the central < 2" where promi- 
nent dust features are present. Moreover, the amplitude 
of the slope is similar to the linear relation of SV'a ver- 
sus Ini? in the bottom- right panel, but with negative 
instead of positive sign. This is consistent with a loga- 
rithmic spiral with the same pitch angle but orientated 
clockwise instead of anti-clockwise, i.e., leading instead 
of trailing. In principle, both leadi ng and trailin g nu- 
clear spirals can exist, as shown bv IWadal (|1994), but 
in their simulations the smaller leading spiral dissolves 
before the long-lived larger trailing spiral fully emerges. 
Alternatively, the spiral-like contribution to the m = 2 
harmonic term, appearing like a "negative image" , might 
result from an asymmetric dust obscuration, mimicking 
a lopsided distortion. 

The often prominent dust lanes along the leading edges 
of bars in galaxies are associated with sho cks in the gas 
strea ming along the length of the bar (e.g. lAthanassoulal 
|1992|) . which in turn lead to velocity j umps across the 

19991 ). Moreover, 
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dust lanes (e.g. iMundell fc Shonel 



meric al models of dust lanes (e.g. 



Gerssen fc Debattistal 



20071) as well as ana,lytica l models of diffuse disks (e.g. 



Valotto fc Giovanellill2004l ) show that dust extinction can 



have a significant effect on the velocity along the line- 
of-sight. Henceforth, we expect the dust and possible 
shocks associated with the spiral features to distort the 
velocity field, but the modeling required to understand 
the specific effects on the non-circular motion is beyond 
the scope of this paper. 

4.3. Feeding the central black hole? 



MEdd = 2.2MQyr- 



Af, 



BH 



VlO^Me 



(23) 



adopting e = 0.1 for the radiative efficiency, and 
a mass A/bh = 1.2 x 10^ M© for the central black 
hole in NGC 1097 — based on the measured central 
stellar velocity dis persion of o-^, = 196 ± 5kms~^ 
(|Lewis fc Eracleoui i2006i ) substituted in the Mbh — 
CT^ relation (|Tremaine et al.l I2002D — yields M^Edd — 
2.7Mq yr~^. This value is indicated in the bottom-right 
panel of Figure [5] by the upper dashed horizontal line. 
The lower dotted horizontal line at M = O.OlMsdd is 
the approximate transition from a Seyfer tl to a L INER 
active galactic nucleus (see for a review lHol[2005l) . and 
thought to be the rate below which mass accretion be- 
comes radiB^tivelj;^^ (see for reviews QuataertI 
[200lt INaravanI [20051) . Whereas NGC 1097 is typically 
classified as a LINER galaxy, monitoring of the nucleus 
reveals evolution in its activity up int o the Seyfertl 
regime ([Storchi-Bergmann et al.ll2003( ). iNemmen et al.l 
find that the observed optical to X-ray spectral 



energy distribution in the nucleus of NGC 1097 is consis- 
tent with an inner radiatively inefficient accretion flow 
plus outer standard thin disk, with a mass accretion rate 
of M ~ 6.4 X 10"'^ A^Edd- These observational constraints 
are fully consistent with our estimate of the mass inflow 
rate, down to il ~ 0.011 Mq yr-i(~ 4.2 x lO^^AfEdd) 
at a distance R — 70 pc (~ 1") from the center, where 
the gas kinematics are still accurately measured and well 
described by the two-arm spiral perturbation model. 

Our mass inflow rate is significantly lower t han M 
O.6M0yr"^ estimated bv IStorchi-Bergmani] (|2007f ) at 
a distance R — 100 pc. In the latter estimate, a non- 
circular motion of 50kms^^ is adopted for the infiow 
velocity, whereas we find from our spiral model that 
even though the contribution of the radial inflow dom- 
inates over elliptic streaming, it is still only a fraction 
(~ 13kms~^ at 100 pc) of the non-circular motions. 
Moreover, the estimated electron density of ~ 500 cm"'^ 
is similar to the average density we derive directly from 
emission line ratios, but more than an order of magni- 
tude higher than the expected (over)density in the spi- 
ral arms (~ 43cm~^ at 100 pc). Finally, instead of 
assuming a fixed opening angle of 20°, we incorporate 
the geometry from the spiral model, and allow for the 
scale height to vary, in order for the disk to remain 
marginally stable. In this way, we obtain an estimate of 
the mass inflow rate fully consistent with the expected 
mass accretion rate onto the central BH in NGC 1097, 
without having to invoke an unclear flUing factor (e.g. 
IStorchi-Bergmann et al.|[2007l ). Even so, the mass inflow 
rate might be even further constrained, in particular by 
measuring emission line ratios at higher signal-to-noise 
to obtain a more accurate estimate of the (over) density, 
as well as using other emission lines to break the degen- 
eracy with temperature. 

Even though the ratios of emission lines constrain the 
total gas density in the nuclear spiral, most likely only 
part of the gas is ionized. On the other hand, the high 
amount of radial damping needed for the spiral model to 
match the observed non-circular motions (Section 13. 2p . 
indicates strong dissipation, possibly due to shocks that 
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in turn can ionize a significant fraction of the gas. At 
the same time, it imphes that this "hot" gas does not 
necessarily show the same kinematics as the underlying 
"cold" gas, as traced for example by Hi and CO measure- 
ments at rad io wavelengths, like those presented recently 
iDavies et alJ ([2009.) . In the latter study, the projected 
molecular gas distribution in the center of NGC 1097 
shows after subtracting elliptical isophotes, a clear two- 
arm spiral structure (their Fig. 3), consistent with a 
m = 2 spiral perturbation in the gravitational potential. 
However, the corresponding velocity field after subtrac- 
tion of an axisymmetric disk model (their Fig. 5) does 
not reveal a clear three-arm spiral structure, as seen in 
the non-circular motions of the ionized gas (Figure H]) . 
It is possible that the kinematics of the ionized gas bet- 
ter trace the nuclear spiral, whereas the corresponding 
three-arm spiral is weaker in the non-circular motions of 
the cold gas, and perhaps more contaminated by a m = 2 
contribution through asymmetric dust obscuration. 

5. SUMMARY AND CONCLUSIONS 

We presented harmonic expansion of the line-of-sight 
velocity field as a suitable method to identify and quan- 
tify possible structures in the non-circular motions, in- 
cluding nuclear spirals. We showed that a weak pertur- 
bation in the gravitational potential of harmonic number 
m, causes the surface brightness to also exhibit m dis- 
tortion, but leads to m — 1 and m + 1 harmonic terms 
in the velocity field. In the case of a m-arm spiral per- 
turbation in the gravitational potential, we found that 
the corresponding (m -I- l)-arm and (m — l)-arm spi- 
rals in the velocity field are respectively less and more 
tightly wound, with pitch angles approximately related 
as m cot Cm ~ (m ± 1) cotCmii- In Appendix |A| we de- 
rived an analytic perturbation model, which allows for 
a simple estimate of the fraction of the measured non- 
circular motions that is due to radial fiow. 

We applied this method to the emission-line veloc- 
ity field within the circumnuclear starforming ring of 
NGC 1097, obtained with the GMOS-IFU spectrograph. 
The resulting non-circular motions reveal a three-arm 
spiral structure, and from the corresponding m = 3 har- 
monic terms, we measure a pitch angle (^3 = 63 ± 3°. 
We linked the m = 1 and m — 3 harmonic terms in 
the velocity field to a weak perturbation of the gravita- 
tional potential due to a two-arm nuclear spiral with an 
inferred pitch angle C2 — 52 ± 4°. We showed that this 
agrees with a two-arm dust spiral in (optical) images of 
the center of NGC 1097, although additional spiral struc- 
ture might be present as a result of higher-order even 
harmonic terms. We argued that it is unlikely that the 
presence of also m = 2 harmonic terms in the velocity 
field is due to lopsidedness in the gravitational poten- 
tial nor in the nuclear spiral itself, which we expect to 



be due to long-lived density waves in the gas driven by 
the bi-symmetric large-scale bar. Instead, we postulated 
that an asymmetric dust obscuration mimics a lopsided 
distortion, and gives rise to the additional even harmonic 
terms in the velocity field. 

To match the measured odd harmonic terms in the 
velocity field of NGC 1097, a spiral perturbation model 
with a large amount of radial damping is required. This 
indicates both strong dissipation and, that, in addition 
to elliptic streaming, a significant fraction of the non- 
circular motions is due to radial inflow. We combined 
the inferred radial inflow velocity with the gas density 
in, and the geometry of, the spiral arms to estimate the 
mass inflow rate as function of radius. Wc calculated 
the gas density from the variation in the flux ratio of 
the [Sii] AA6716.4,6730.8 emission-line doublet, and we 
incorporated in the geometry the scale-height assuming a 
marginally stable disk. The resulting mass inflow rate de- 
creases to M ~ 0.011 Mq yr~^ at a distance of R — 70 pc 
from the center. We showed that this corresponds to 
M ~ 4.2 X 10~^ MEdd, where the latter Eddington accre- 
tion rate is onto a black hole in NGC 1097 with a mass 
Mbh = 1-2 X 10^ M0 based on the central stellar veloc- 
ity dispersion. This rate is fully consistent in general 
with the active galactic nucleus in NGC 1097 varying be- 
tween LINER and Seyfertl, and in particular with previ- 
ous mass accretion models fitted to the observed spectral 
energy distribution in its nucleus. 

We conclude that the line-of-sight velocity can provide 
not only a cleaner view on nuclear spirals than does the 
associated dust, but that the presented method also al- 
lows one to quantitatively study these possibly important 
links in fueling the centers of galaxies, including a handle 
on the mass infiow rate as a function of radius. 
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APPENDIX 

GASEOUS ORBITS IN A WEAKLY PERTURBED GRAVITATIONAL POTENTIAL 

We derive analytic solutions of gaseous orbits in a weakly perturbed gravitat ional potential (for modelin g of a strong 
perturbation see e.g. ISpekkens fc Sellwood 2007). Wc follow the tr eatment of iBinnev fc Tremaind (Il987l . p. 146-148) 
for a weak harmonic perturbation. Like lSchoenmakers et al.l ()1997f ). we include an additional phase-shift 4>m{R) as in 
equation ([5]) to accommodate for a spiral perturbation. Since gas at the same spatial location has the same velocity, 
we look for closed loop orbits that do not intersect themselves. However, to take into account the dissipative nature 




Fig. 5. — An axisymmetric logarithmic potential perturbed by a weak m = 2 logarithmic spiral with pitch angle (^2 = 50°. The top 
panels show for an inclination i = 30°, projected maps on the sky-plane {x',y') of the line-of-sight velocity, non-circular motions, and 
separately, the contribution of the m — 1 and m + I harmonic terms. The first panels in the middle and bottom rows show the intrinsic 
maps in the equatorial plane (x, y) of respectively the azimuthal and radial velocity due to the spiral perturbation. The overplotted open 
circles and stars trace the maxima as function of radius, resulting in the radial profiles shown in the second panel in the bottom row. The 
second panel in the middle row shows the radial profiles of the circular velocity Vc of the axisymmetric potential, while the dotted curve 
shows ci, which includes the contribution from the spiral perturbation. Moreover, the falling dashed curve is the angular frequency H, 
which intersects the horizontal dashed line of the assumed pattern speed Qp, at the corotation radius Re- The latter is indicated in all 
panels either by a vertical dashed line or dashed ring on the maps. The third and fourth panels in the middle row show the radial profiles of 
the harmonic coefficients c^ii and Sm±i, divided by sini to correct for inclination. The third and fourth panels in the bottom row show 
the corresponding phase shifts 4'm±i, defined in equation (|4]l, against Ini?. In this way, a spiral structure shows up as a smooth variation, 
which becomes nearly linear because the assumed spiral perturbation is logarithmic as in equation l|12|l . The slopes of the fitted solid lines 
yield the indicated values of the pitch angles Cm±l of the (m_± l)-arm spirals shown in the corresponding projected maps in the third and 
fourth panels in the top row. In turn, based on equation (1151 . both values provide an estimate as indicated of the pitch angle £^2 of the 
m = 2 spiral perturbation. The mild difference with the input value C2 = 50° is due to the radial variation of the additional terms 9m±i 
indicated by the dot-dashed lines. 



of gas, we include, like I Wadal (|1994r ). radial damping, so that the first-order equation of motion becomes 
Ri + 2XkoRi + kIRi = -Rq (A cos 7/ + Bsinri) , rj ^ m[0o(i) - '/'m(-Ro)], 



(Al) 



where A controls the amount of radial damping. The subscript zero refers to zeroth-order with constant i?o and 
00 (i) — {^Q — ^p)t, with Hp the pattern speed of the perturbation. We have introduced the epicycle frequency, 
= 2f2^(l -|- a), and angular frequency, fl = v^/R, with a = dlnuc/dlni? the logarithmic slope of the circular 
velocity, v'^{R) — d$o(-R)/dlni?, of the axisymmetric part of the gravitational potential. Furthermore, 



i?2 



dln$^ 



1-UJ 



dlnR 



B = cot C„ 



(A2) 
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where we have introduced uj — flp/il, and is the pitch angle as defined in equation (jlip . Solving equation (lAip we 
find that the solution for closed loop orbits is 

i? = i?o[l ^ (acos77 + 6sin?7)], (A3) 

(j) = (j)o + [2(asin77 — bcosrj) — ^sin77]/[m(l — uj)] (A4) 

VB. = Vc {m{\ — ix')[osinr7 — fecosry]} , (A5) 

= fc {1 + (1 + a) [(a - cos 77 + 6 sin 77)} , (A6) 

with a appearing in the last line because of a first-order conversion from guiding center [Ro,(f>o) to a point (i?, 0) in 
the observed velocity field. Moreover, 

AA-BA BA + AA . /Ay^ 

where we have defined A — k'^ — m^{il — ilp)^ and A — 2\Km{Vl — VLp). The orbit solutions have a singularity at the 
corotation radius where SI = f7p (w = 1), because the adopted epicycle approximation breaks down. Without radial 
damping (A = A = 0), the (coUisionless) orbit solutions also have singularities when A = 0, i.e., at the Lindblad 
resonances given by — n/m = Qp. 

Next, we assume to first order (j)Q ^ (j) and replace it hy ip = (j) — (pios + 7r/2 which is zero along the line of nodes. We 
also define tpm = 4>m — <P\os + 7''/2, so that 77 — mlip — ipm{R)]- This allows us to recast the above expressions for vr 
and in multiple angles of "0 as in equations ([7|) and ([8]), with 

Ci^ = — 771(1 — uj)[a sm{m(p,n) + bcos{mtpjn)], (A8) 

Sfl^77i(l — uj)[a cos{m(pm) — bsm{m(pm)], (^9) 

C0 = (1 + a){a - ^) cos{mipm) - (1 + a)b sm{mipm) ■ (AlO) 

s^ — {l + a){a — £_) sin{mipm) + (1 + a)b cos{mipjn) , (All) 

Substitution in equation ^ results in the expression for the line-of-sight velocity in equation ([5]), where the harmonic 
coefficients are given by 

c™±i=A±cos(77zi^„) - B±sin(77i(^„) (A12) 
Sm±l = A± Sin(771(/?m) + B± cos(r7i(/?m), (A13) 

with 

A± = ^V. {[{l + a)Tmil-u;)]a-{l + a)^}, B± = ^¥4(1 + a) T mil ^ Lj)]b. (A14) 

Substituting A± — Km±i cos0m±i and B± = Km±i sin0m±i, with 

i^^±i=A|+S2, tan^^ii =B±M±, (A15) 

we arrive at the form given in equation p3p . 

In Figure[5l we present an example of a weakly perturbed axisymmetric logarithmic potential defined in equation (|19p , 
with vq = 250 kms""'^, Rc = l.Okpc and /3 = 0. The amplitude of the perturbation is a factor Cp = 0.02 times the 
axisymmetric logarithmic potential, while the angular dependence is due to a logarithmic spiral defined in equation (fT2l) , 
with pitch angle C2 = 50°. This results in non-circular motions which contribute a 1-arm and 3-arm spiral structure 
to the observed line-of-sight velocity field as shown in Figure [5] for an adopted inclination of i = 30°. 

In general, Km±i and 0m±i in equation (|A15p above depend in a rat her complicated way on the gr avitational 
potential. However, in case of the linear spiral density-wave theory (|Shu et al.l U973I: iCanzian fc Allenl Il997.) the 
expressions reduce to 



Km±i = 



K 777(ri — flp) 

2n ^ K 



tan6l„±i = ± cot Cm, (AI6) 



where the constant amplitude Vm measures the strength of a tightly wound spiral without radial damping. In this 
case, cot Cm ^ 1 and A = A = 0, so that 6 = B/A ^ a = A/ A, and b ^ in equation (|A7p . and hence tan6'm±i = 
B/A cx cot Cm and Km±i — B±. The expression for B± in equation (|A14p is proportional to the reduced expression 
for Km±i in equation (|A16p . since by substituting k'^ — 2fl^{l + a) and lu — ilp/il, the term in square brackets in 
equation (|A16p is proportional to [(1 + a) ^ m{l — uj)] in equation (jA14p . 
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